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ABSTRACT 

We present new calculations of X-ray polarization from accreting black holes 
(BHs), using a Monte-Carlo ray-tracing code in full general relativity. In our 
model, an optically thick disk in the BH equatorial plane produces thermal seed 
photons with polarization oriented parallel to the disk surface. These seed pho- 
tons are then inverse-Compton scattered through a hot (but thermal) corona, 
producing a hard X-ray power-law spectrum. We consider three different mod- 
els for the corona geometry: a wedge "sandwich" with aspect ratio H/R and 
vertically-integrated optical depth tq constant throughout the disk; an inhomo- 
geneous "clumpy" corona with a finite number of hot clouds distributed randomly 
above the disk within a wedge geometry; and a spherical corona of uniform den- 
sity, centered on the BH and surrounded by a truncated thermal disk with inner 
radius Redge- In all cases we find a characteristic transition from horizontal po- 
larization at low energies to vertical polarization above the thermal peak; the 
vertical direction is defined as the projection of the BH spin axis on the plane of 
the sky. We show how the details of the spectropolarization signal can be used 
to distinguish between these models and infer various properties of the corona 
and BH. Although the bulk of this paper focuses on stellar-mass BHs, we also 
consider the effects of coronal scattering on the X-ray polarization signal from 
supermassive BHs in active galactic nuclei. 



Subject headings: black hole physics - accretion disks - X-rays:binaries 



- 2 - 



INTRODUCTION 



A recent flurry of new mission proposals has renewed interest in measuring X-ray polar- 
ization from a variety of astrophysical sources. The Gravity and Extreme Magnetism SMEX 
( GEMS) missioiJl], which has recently been approved for funding in the latest round of NASA 
Small Explorer prop osals, should be able to detect a degree of polarization S < 1% for a 
flux of a few mCrab (IBlack et al.ll2003l : iBellazzini et al.ll2006l : ISwank et al.ll2009l ). A similar 
detector for the International X-ray Observatory (IXO) could achieve sensitivity roughly 
lOx greater (< 0.1% degree of polarization; Jahoda et al. 2007, Costa et al. 2008). Projects 
like these could potentially detect a large number of galactic and extra-galactic sources at 
the 5 ~ 1% level, including stellar- mass black holes, magnetars, pulsar wind nebulae, and 
active galactic nuclei. In this paper, we focus on accreting b lack holes (BHs) in the "Hard" 
or "Steep Power Law" state (IRemillard fc McClintockll2006l ). which are characterized by a 
broad-band spectrum with a thermal peak around 1 keV and a strong high-energy power- 
law component extending above 100 keV. Early estimates suggest that the typical level of 
polarization from these sources should be a few percent in the 1 — 10 keV ra nge, depending 



on the geometry of the accretion system and the inclination to the observer (jConnors et al. 



1980 : Sunyaev &: TitarchuklEoSsI ). We also extend these results to AGN, whose hard X-ray 
properties are in many ways qualitatively similar. 

Symmetry demands that in the Newtonian limit the observed polarization from a fiat 
disk must be either parallel or perpendicular to its rotation axis. However, the effects of 
relativistic beaming, gravitational lensing, and gravito-magnetic frame- dragging can break 
that symmetry and give a non-trivial net rotation to the integrated polarization vector. Be- 
cause the temperature in an accretion disk should increase closer to the BH, where these 
relativistic effects are strongest, it was predicted long ago that the observe d angle and degree 
of polarization of therrnal disk emission should d epend on photon energy (jStark fc Connors 
19771 : IConnors fc Starklll977l : IConnors et al.lll980l ). In more recent years, discussion of polar- 
ization in accreting black holes has expanded to include a numb er of other aspects of thermal 



disk emission, such as UV and X-ray emission from AGN disks (ILaor et al.lll990l : iMatt et al. 



19931 ) and "lamp post" models for ir radiating the acc retion disk wit h a no n-thermal source 
above the plane flDovciak et al.ll2004l ). Quite recently, iDovciak et al.l (120081) investigated the 
effect of atmospheric optical depth on the disk's polarization signal, and lLi et al.l (120081 ) ap- 
plied the original calculations of therm al X-ray polarizatio n to the problem of measuring the 
inclination of the inner accretion disk. iDavis et al.l (120091 ). using a Monte Carlo ray-tracing 
code on data from shearing-box simulations, estimated the effects of magnetic fields on the 
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polarization of thermal emission. ISchnittman fc KrolikI fl2009[) fhereafter SK09) showed the 
importance of including returning radiation (j Cunningham! 1 1 9 76l ) when calculating the p olar- 
ization of BHs in the thermal state, an effect previously noted by lAgol fc KrolikI (120001 ) . 



In addition to the relatively well-understood thermal state, it has also long been known 
that most AGN and stellar-mass BHs can produce significant levels of harder X-rays with 
energies well above the thermal peak. In both cases, the hard flux is thought to be produced 
via i nverse Compton scat t ering of the disk photons in a corona of hot (yet thermal) electrons 



(e.g. iHaardt fc Maraschil (119931 )). From the shape of t his hard spectrum, the basic geometry 
and optical depth of the corona may be con strained (IHaardt et al.l Il994l : iPietrini Sz Krolik 
I995I : Istern et al.l[l995l : IPoutanen et aDll997h . 



Despite these constraints, there remain sizable uncertainties about the nature and ge- 
ometry of the hard X-ray emitting region. Popular models for the hard state of stellar-mass 
BHs include a cool disk truncated at large radius ('~ lOOM) surrounding a hot, radiatively 
ineffi cient flow Jcierlihski et al.lll997l : [McClintock et al.lboOllJEsin et al.ll200ll : lDone fc Trigo 
2OO9I ). or alternatively a more ex tended disk surrounded by an optically thick hot corona. 



possi bly in the form of a hot wind (IBlandford fc Begelmanll2004l : iMiller et al.ll2006l : iReis et al. 



2OO9I ). Even less is known about the steep power-law (SPL) state. As for the hard state, 
most popular models are based on the invers e Compton scattering of ther mal seed photons 
from a thin disk surrounded by a hot corona (IZdziarski Sz Gierlinskil 120041). but bulk Comp- 



tonization of a cony erging accretion flow has also been suggested (ITitarchuk fc Shraderll2002 



TuroUa et al]l2002h . Although the geometry of the corona could be as simple as a uniform 
slab in the SPL state (IZdziarski et al.ll2005l ). the hard state of galactic BH binaries, as well 
as the X-ray emission from AGN, ar e more likely caused by clumpy, inhomogeneou s coronae, 
possibly caused by magnetic flares (IHaardt et al.lll994l : iPoutanen fc Fabianlll999l ). 



Here we explore the X-ray polarization signatures of three simple models for the corona 
geometry: a smooth sandwich with uniform optical depth in the vertical direction; an inho- 
mogeneous model made of a finite number of spherical clouds, randomly distributed above 
the disk; and a truncated thin disk surrounding a spherical corona. In all cases, seed photons 
are emitted from the thermal, optically thick disk and are up-scattered in the hot corona (in 
the truncated disk case, we also include low-energy seed photons embedded in the corona). 
We find that in all cases, the spectrum comprises a thermal peak (~ 1 — 3 keV for stellar- 
mass BHs; ~ 30 — 100 eV for AGN) and a power-law component dominating above that, 
which typically makes up ~ 50 — 80% of the total fiux. The gross features of the polarization 
spectra are quite robust: horizontal polarization with a few percent amplitude at low ener- 
gies, with a transition to vertical orientation above the thermal peak, where the polarization 
amplitude can be as large as ~ 10%. The specific details of the polarization spectrum, i.e.. 



-4- 



the amplitude of polarization at low and high energies, and the shape and location of the 
transition, provide constraints on the global geometry, temperature, and optical depth of the 
corona, as well as the BH mass, spin, accretion rate, and the observer inclination angle. 

In this paper, we begin in Section [2] with a brief overview of the computational methods 
used in the calculations. In Sections [3H5] we present our results for three different corona 
geometries: sandwich, hot spots/clumps, and a sphere embedded in a truncated disk, all in 
the context of galactic X-ray binaries. In Sectional we apply the results to AGN. In Section 
[7] we discuss possible applications to observations. 



2. METHODOLOGY 



Al though a detailed descripti on of the ray-tracing code will be given in a companion 
paper (jSchnittman fc KrolikI |2010| ). we give a brief summary of the relevant physics here. 
The basic geometry is described by a thin disk with inner radius Redge, typically placed at 
the inner-most stable circular orbit (ISCO), but possibly at smaller radii so long as -Redge is 
outside the horizon. In the latter case, the gas inside the ISCO follows plunging trajectories 
along geodesies with energy and angular momentum determined by the ISCO values. The 
orbital angular momentum of the disk is prograde and aligned with the BH spin axis. Seed 
photons are emitted from the t hin disk with a diluted bla ck-body spectrum characterized 
by a hardening factor / = 1.8 (IShimura &: Takaharal 19951) . The local flux is given by the 
Novikov-Thorne emissivity profile (INovikov &: Thornelll973l ). but can be modified to include 
emission inside the ISCO, as described in SK09. 

The seed photons have initial polarization parallel to the disk surface in the local 
f luid frame, as determ ined by the classical result for scattering-dominated atmospheres 
( 1 Chandr asekharl 119601 ) . The degree of polarization varies from zero for photons emitted 
normal to the disk surface up to ~ 12% for an inclination angle of 90°. In addition to the 
polarization effects, the scattering of the outgoing flux causes limb-darkening, effectively 
focusing the emitted radiation in the direction normal to the disk surface. We define our 
seed photons in accordance with the polariz ation and limb-darkeni ng factors tabulated as a 
function of emission angle in Table XXIV in I Chandr asekharl (jl960l ) . 



After leaving the disk, the photon packets follow geodesic trajectories around the black 
hole, eventually reaching a distant observer, returning to the disk via gravitational deflection, 
or getting captured by the horizon. En route, many photons scatter in the corona. We 
describe these events with the classical Thomson electron cross section. The differential 
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probability for scattering along a path length dl is simply 



scat 



1 -e- 



-dr 



1-e 



-Kp dl 



where k = 0.4 cm^ g~^ is the opacity to electron scattering, p is the local mass density of 
the corona, and r is the optical depth . Unl ike some recent Monte Carlo ray-tracing codes 
(e.g. iDavis et al.l (120091 ) : lDolence et al.l (120091 )) that use the more physically accurate, energy- 



dependent Klein-Nishina cross section, the Thomson cross section allow s for the use of photon 
packet s that include the entire broad-band spectrum, as described in ISchnittman fc Krolik 
( 120 lOl ). This simplification gives improved computational efficiency, albeit at the loss of 
physical accuracy for above ~ 200 keV, where the classical and relativistic electron cross 
sections begin to diverge. For most stellar-mass BH sources, and for most X-ray polarimetry 
missions in the foreseeable future, it should be quite safe to focus on photons with energies 
less than 100 keV. 

At each step along the photon geodesic, the probability for scattering is calculated 
according to equation ([T]), and then a uniform random number in [0, 1) determines whether 
that photon scatters at that location!^. When the photon does scatter, we first transform 
the photon 4-momentum k and polarization vector f into the local inertial frame of the 
corona, then do a special-relativistic boost into the electron's rest frame, where the electron 
velocity is taken from an isotropic thermal distribution. In the electron's rest frame, the 
scattering event is treated completely classically, with a differential cross section given by 
jRvbicki fc Lightman|[l979l ) 



dn ~ 2^^°^^^ 



5) (cos^ e + l) + 2S cos^ e cos^ tp + 25 sin^ 



(2) 



where tq = 2.82 x 10"^'^ cm is the classical electron radius, 6 is the degree of polarization, 
is the angle between incident and outgoing photon directions, and ip is the polarization 
angle, measured with respect to the scattering plane. The angles 6 and ip are selected 
appropriately from the probability distribution corresponding to ([2]). The outgoing Stokes 
parameters are given by 



Q' 
u' 



-(/iicos^e 



cos^ Q-h 
U cos 9 , 



(3a) 
(3b) 
(3c) 



^Our adaptive-step Cash-Karp geodesic integrator takes smaller steps where the electron density is high, 
ensuring that dr « 1. 
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where /y and I± are the components of the intensity with polarization parallel to and per- 
pendicular to the scattering plane, respectively. These Stokes parameters are then used to 
reconstruct the new polarization degree and angle through 

X = Q/I, (4a) 
Y = U/I, (4b) 



and 



6 = {X^ + yY', (5a) 
= -tan^\Y/X). (5b) 



The scattered photon packet is boosted back to the coronal fluid frame, then transformed 
back into the coordinate basis and continues propagating along its new geodesic. This 
"transform-boost-scatter-boost-transform" algorithm automatically carries out the inverse 
Compton change in photon energy, giving the outgoing photon an average energy increase 
of 7^, where 7 is the electron Lorentz factor (7 1 + kT/{meC^) for thermal electrons). 
The polarization 4-vector is parallel-transported simply by satisfying the constraints k • 
f = and f ■ f = 1 and using the complex- valued Walker-Penrose integr al of motion k-wp 



(I Walker fc Penrose! 1 19701 ). Analogous to the way that Carter's constant (jCarterlll968l ) can 
be used to constrain 4-velocity components, kwp can be used to reconstruct the polarization 
vector at any point along the geodesic. The two orthonormality conditions stated above, 
along with the real and imaginary parts of kwp, give a total of four equations for the four 
components of the polarization vector f . 

Some photons encounter the disk before reaching their final destination. When we treat 
models relevant to AGN, photons striking the disk are absorbed because our primary interest 
is in photons < 10 keV, and AGN disks are generally thought to have enough photo-ionization 
opacity to have little albedo at these energies. In the stellar-mass black hole case, disks are 
expected to be h ighly reflective. For these models, we use the results given in Section 70.3 of 



Chandrasekharl (jl960l ) for diffuse reflection, i.e. multiple scattering events in a semi-infinite 



plane. 

When a ray reaches an observer at infinity, the polarization vector is projected onto 
the detector plane. In our convention, the y-axis of the detector is parallel to the projected 
symmetry axis of the black hole and accretion disk, so ^/^ = 0° corresponds to "horizontal" 
polarization and if) = ±90° is "vertical." Finally, we integrate over all photon bundles to 
obtain the energy-dependent Stokes parameters 1^, Q^, and f/^, which in turn give 6u and 
ipi,, our preferred observables for most of the results presented below. 
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3. WEDGE GEOMETRY 

We begin by considering the simplest corona geometry, an isothermal layer of uniform 
vertical optical depth forming a wedge with constant opening angle tan 9c = H/ R. At each 
point in the disk, the atmosphere has an exponential profile in the vertical direction with 
p{z,R) = poexp[—z/H{R)]. The vertically-integrated optical depth tq is constant, giving 

We set the local rest frame of the corona to be corotating with the underlying disk: = 
= 0, u'^ /v} = fl{R), and u^u^ = —1. Q{R) is the orbital frequency (measured at infinity) 
of a planar circular orbit at radius R, which in Boyer-Lindquist coordinates is 

^^^^ = {R/Mrl + a/M ■ 

Figured] gives a schematic view of the accretion geometry, showing how the seed photons 
originate in the midplane, then scatter through the hot corona before reaching a distant 
observer. The coronal scattering has two major effects on the observed polarization signature. 
First, it changes the underlying spectrum by inverse Compton scattering a portion of the 
photons, boosting them to higher energies, and thus producing a harder spectrum. Second, 
it changes the amplitude and orientation of the net polarization at high energies, rotating 
it from horizontal to vertical, similar to the way returning radiation gets scattered into 
a vertical orientation by a thermal disk (S K09). The coronal polarizati on rotation effect 



is essentially the same as that described in ISunyaev fc TitarchukI (Il985l ). who considered 



the up-scattering of low-energy photons embedded in a hot electron corona with planar 
geometry, finding strong vertical polarization for high-energy photons and observers at high 
inclinations. 

The cause of this rotation in systems with moderate optical depth (tq < 2) can be 
understood easily from the scattering geometry in Figure [TJ Photons initially emitted in 
directions near the disk plane are almost certain to scatter because they face a large optical 
depth. To reach an observer who views the disk nearly edge-on, they must stay in the disk 
plane even after scattering, but this requires a vertical polarization direction. Additionally, 
those photons that scatter multiple times in the corona — and are thus boosted to higher 
energies — are geometrically more likely to move in the plane of the corona, parallel to the 
disk. This further increases the amplitude of their vertical polarization and preferentially 
scatters them to infinity with large emission angle, leading to a limb- brightening effect at high 



energies with respect to the classical Chandrasekhar result (ISunyaev fc TitarchukI 119851 ). 
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For hot coronae with high optical depth, the multiple scatters cause the polarization 
and limb-darkening to tend towards the Chandrasekhar limit, but with the photon energy 
dependence given by a Wien spectrum. In the limit of a cold corona, where there is no 
transfer of energy via inverse Compton scattering, there is no way to distinguish between 
photons coming directly from the disk atmosphere and those that are scattered in the corona, 
and the classical result is reproduced regardless of the coronal optical depth. 

To improve our qualitative understanding of the polarization effects of the global corona 
geometry, we plot in Figure [2] a series of images of the BH accretion disk plus corona, sorted 
by the scattering history of each photon. Figure [2^ shows only photons that travel directly 
from the disk to the observer; Figure [2]d shows flux from return radiation — photons that are 
bent by gravitational lensing and scatter off the disk to the observer; Figure [2]3 is the flux from 
photons that scatter exactly once in the corona; and Figure is made of photons that scatter 
multiple times in the corona. The observed intensity is color-coded on a logarithmic scale 
(normalized to the peak intensity from all photons) and the energy-integrated polarization 
vectors are projected onto the image plane with lengths proportional to the local degree of 
polarization. The black hole has a spin of a/M = 0.9, mass M = lOM©, and accretion rate 
such that the thermal flux alone totals 10% of the Eddington luminosity. The corona has a 
temperature of 100 keV, scale height H/R = 0.1, and vertical optical depth tq = 1. These 
parameters lead to a coronal luminosity 2 — 3 times that of the thermal luminosity. The disk 
is rotating in the counter-clockwise direction, and the observer is located at an inclination 
of 75° to the rotation axis. The intensity maximum at the left of each image is caused by 
relativistic beaming and Doppler boosting. Photons from the far side of the disk are bent 
by gravitational lensing, giving the appearance of a warped disk with the far side having a 
smaller effective inclination angle. 

Not surprisingly, the direct image looks quite similar to that for the direct radiation from 
a simple thermal disk with no corona (SK09). The polarization is modest and predominantly 
horizontal, with relativistic effects most noticeable in the inner disk, where beaming and 
lensing combine to rotate the observed angle of polarization at higher energies. The degree 
of polarization is reduced above and to the left of center, where the effective inclination of 
the disk is smaller, while the opposite occurs in the receding section of the disk on the right 
side of the image. However, when carefully compared to the non-corona system (Fig. 1 of 
SK09), we see a clockwise shift in the peak flux distribution because the level of direct flux 
from a given patch of the disk is dependent on the optical depth to scattering along the 
geodesic connecting that patch with the observer. For the wedge geometry, this means that 
where the effective inclination is smallest, the observed flux is greatest, leading to additional 
limb-darkening in the direct component. 
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The return radiation image is also similar to that for the thermal disk alone, although 
now the presence of a scattering corona strongly suppresses any return flux from large radii. 
Since most returning photons pass very close to the light orbit around R = 3M, they would 
necessarily intersect the outer disk at very high angles of incidence, and thus face very large 
optical depth through the wedge corona. As shown in SK09, the flux from returning ra- 
diation, while small compared to the total flux, is very highly polarized, and in a vertical 
direction, perpendicular to that of the direct flux. Since the majority of the returning radia- 
tion originates in the hot inner regions of the disk, there is a dependence of the polarization 
direction on photon energy. As can be seen from Figure [2)d, the angle of polarization for the 
return radiation is negative, as measured in the image plane, with ip = defined along the 
X-axis. This can also be understood from the global geometry of the system, as shown in 
Figure [TJ Photons emitted from the right side of the disk, deflected by the BH and incident 
on the left side at some finite angle, will scatter towards a high-inclination observer with 
large polarization in the ip < orientation. Of course, the opposite occurs for photons emit- 
ted from the left side and scattered on the right side, but they get preferentially scattered 
away from the observer by relativistic beaming. 

In Figure ^ we show the flux from photons experiencing a single coronal scatter. Since 
every photon in this image scatters exactly once, there is no preference for oblique emission, 
and photons emitted at all inclination angles are equally weighted. If anything, for tq = 1, 
the seed photons emitted normal to the disk surface are more likely to scatter exactly once 
and then escape to an observer at high inclination, thus explaining the prevailing horizontal 
polarization seen here. Lastly, Figure shows the multiply-scattered photons. These have 
strong net vertical polarization due to the geometric effects described above and shown 
schematically in Figure [H 

Unlike the returning radiation, the coronal flux has a net positive polarization angle 
ip > when viewed edge-on. This can again be understood from the schematic image of 
the wedge corona in Figure [H It is clear from this geometry that for two photons, both 
emitted in the plane of the sky at the same angle icm to the disk normal, the photon emitted 
in the direction away from the black hole will have a greater path length within the wedge 
corona. However, due to the increasing scale height and constant optical depth in the vertical 
direction, it will pass through lower-density gas. Combining these two effects — path length 
and coronal density as a function of emission angle — we find that the optical depths for 
ingoing and outgoing photons are 
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and 

rout(^cm) = — p + — • 8b) 

2 LcOS(Zem + yc) COS 2cm. 

Here 6c is the opening angle of the corona and we restrict Zgm < 7r/2 — so that the outgoing 
optical depth is finite. 

Since Tout > Tin for each Zem, the "tilt" of the polarization vector for corona scattering 
is opposite that of return radiation. For the region of the disk beamed towards the observer 
(left side of images in Fig. [2]), this means that photons emitted in the direction away from 
the BH are more likely to be scattered in the wedge corona, giving a net polarization oriented 
with ip > 0, as can be seen in the lower panels of Figure [21 In practice, we can't ever really 
know whether a given photon came directly from the disk or was scattered in the corona or 
off the disk, but the images do provide important qualitative understanding of the geometric 
effects involved. 

For a more quantitative picture, we plot in Figure [3] the observed flux, polarization 
degree, and polarization angle as a function of energy for observer inclinations of i = 45°, 60°, 
and 75°. We sort the photons into direct (dotted curves), scattered (dot-dashed curves), and 
total flux (solid curves). For a BH mass of IOMq accreting at a rate ~ O.lLgdd; the thermal 
peak is near 1 keV, with essentially no flux from the disk above ~ 10 keV. With tq = 1 and 
Tc = 100 keV, the Compton y pararneter i s of order unity, so that the up-scattered spectrum 



has -Fjy oc u'^ (IRybicki fc Lightmarull979l ). As the inclination of the observer increases, the 
optical depth to the disk increases, thus lowering the relative contribution from the direct 
radiation. In all cases, the up-scattered radiation dominates the flux above a few keV. 



The direct radiation behaves quite similarly to that described in IConnors et al.l (jl980l ). 
reproducing the Chandrasekhar result at low energy, and then decreasing in fractional po- 
larization at higher energies as the relativistic effects of the inner disk begin to dominate. 
As discussed in SK09, the angular rotation in ip is greatest for small observer inclinations 
simply because there is a smaller degree of "classical" polarization to overcome. 

When considering only the scattered photons, the polarization is still given by the 
Chandrasekhar limit at low energies. The scattered photons below the thermal peak have 
typically scattered only once in the corona, thus preserving their initial horizontal orientation 
(see Fig. At higher energies, we begin to sample the multiply-scattered photons that 
make up Figure [2]i, giving strong vertical polarization. Thus there is a transition from 
horizontal to vertical orientation even when considering the scattered photons alone. When 
including all the flux (solid curves), this transition point is shifted to somewhat higher 
energies than for the scattered photons alone due to the horizontal contribution of the direct 
thermal radiation. 
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As mentioned above, the effects of return radiation and coronal scattering on the sign of 
the polarization angle are opposite: return radiation leads to a negative angle and the coronal 
scattering in a wedge geometry leads to a positive angle in the transition from horizontal to 
vertical polarization. However, some photons that scatter in the corona are also deflected 
by gravitational lensing and then scattered in the corona on the far side of the disk, leading 
to some ambiguity as to whether they should be classified as return photons or coronal 
photons. From Figure [3] we see that the relative contributions from these two effects are 
also dependent on the observer inclination, with low-inclination systems giving a negative 
tilt and high-inclination leading to a positive tilt. This distinction may be observable in 
practice by comparing the polarization in the hard state to that of the soft state, where 
the transition direction is independent of inclination and therefore can be used to define the 
negative orientation (SK09). 

Using spectropolarimetry to infer physical properties of the source is certainly more 
complicated in the hard state than in the thermal state, for the simple reason that there 
are more free parameters in the underlying model, creating degeneracies in fitting the data. 
Instead of attempting to quantify our ability to measure any particular property of the 
BH system, we present here the effects of separately changing individual parameters with 
respect to a fiducial model with M = IOMq, a/M = 0.9, Lthcrm = O.lLEdd, i = 75°, Tq = 1.0, 
Tc = 100 keV, and H/R = 0.1. 

In Figure H] we show the polarization degree and angle for the fiducial model, varying 
the luminosity in the thermal flux component from O.OlLEdd to L-^dd- Since in all cases the 
thermal peak is well below the corona temperature, the transition energy scales linearly with 
the disk temperature, which is in turn proportional to the quarter power of the luminosity. 
With multiple observations of a single source accreting at various rates as in Figure HJ one 
could confirm t hat the thermal part of t he spectrum is indeed coming from a thin disk with 



constant -Redge (IGierlihski &: Donell2005l ). 



Holding the net thermal flux fixed, but varying the spin, changes the emissivity profile of 
the disk, concentrating more of the emission close to the BH for larger spins (and thus smaller 
Risco)- As in the thermal state (SK09), this leads to a greater fraction of return radiation, 
and thus the transition from horizontal to vertical polarization occurs at lower energy, as 
shown in Figure O In the hard/SPL states, the presence of a corona provides another, 
complementary mechanism for effecting the polarization transition above the thermal peak, 
regardless of the spin parameter. Thus the dependence of the polarization signal on BH spin 
is weaker than in the thermal state. In this way, the sandwich corona acts as a veil around 
the inner disk, obstructing our view of the plunging region where the orbital dynamics are 
most sensitive to the spin. 
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In Figure O we simultaneously vary the coronal optical depth and temperature, main- 
taining a roughly constant Compton y-parameter. In the wedge corona model, the transition 
from horizontal to vertical polarization is largely due to local geometric effects (as opposed 
to more global effects like returning radiation). In this case, the optical depth and opening 
angle of the corona largely determine the scattering history of each photon, and thus the 
polarization degree and angle of the outgoing radiation. Holding the size and shape of the 
corona constant, two photon packets experiencing the same number of scattering events will 
emerge on average with the same polarization signature. Thus the energy of the transition 
is primarily a function of the scattering electron temperature, and hotter coronae lead to a 
higher-energy transition. Similarly, to reach a given energy above the thermal peak, photons 
have to scatter more times in a cooler corona, resulting in a more constrained scattering 
history, and thus higher polarization. At the same time, when changing the optical depth 
of the corona, the overall scattering geometry does change slightly, so the shape of the tran- 
sition also varies in Figure [61 Furthermore, for small optical depths, a larger fraction of 
the photons can return to the disk, giving the negative polarization angle characteristic of 
returning radiation. 

Lastly, we show in Figure [7] the effect of varying the scale height of the corona. Here, 
the fact that return radiation and coronal scattering rotate the polarization angle in opposite 
senses away from horizontal leads to a complicated dependence. The geometrically thinnest 
coronae behave similarly to the razor-thin thermal disks described in SK09, giving a more 
gradual transition from horizontal to vertical with ^/^ < (J^l. As the opening angle of the 
corona increases, the transition flips to ^/^ > as the coronal scattering effects of equation 
([8]) becomes more important. For very large H/R, the lower density of the corona allows 
longer path lengths and photons can sample a larger volume of the accretion flow. At low 
energies this loosening of the geometrical constraints leads to a deviation from the classical 
planar scattering atmosphere. At higher energies, it leads to a higher effective temperature 
as the bulk velocity of the corona (which can be a substantial fraction of c) is added to 
the thermal velocities of the scattering electrons. The combination of these various effects 
eliminates the existence of any simple trends in Figure [71 



■^From cqn. ([8|) we see that when 0c —> 0, there is no difference between photons emitted towards or away 
from the BH, leaving only the global effects of returning radiation. 
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INHOMOGENEOUS GEOMETRY 



Motivated by spectral mode ls of AGN Stern et al. ( 19951 ). whose reasoning also applies 
to galactic BHs in the hard state iPoutanen et al.l (jl997l ). we also consider an inhomogeneous 
geometry where the hot coronal plasma is clumped into a large number of small, dense 
regions. These clouds are distributed randomly above the disk, roughly following a wedge 
geometry. For precisely the same reasons that originally motivated the clumpy corona model, 
it is not generally possible to reproduce the same spectrum as in the sandwich geometry. 
Yet we find that by fixing the electron temperature and conserving the total mass in the 
corona, the net flux in hard X-rays (> 10 keV) is roughly the same. 

We describe the inhomogeneous corona with five parameters: aspect ratio H/R, tem- 
perature Tc, mean vertical optical depth tq, number of clouds per unit radius 2nc {uc above 
the plane and ric below the plane; Uc = dNc/dR is constant throughout the disk, giving a 
larger number of clouds per unit area in the inner disk), and an overdensity factor Pc/Po? 
where Pq{R) is the mean density of a wedge corona with the same scale height and optical 
depth. The distribution of clouds is chosen so as to reproduce the wedge corona in the limit 
of Uc ^ oo and Pc/ Po ~^ 1- To do this, we require that the average mass in the inhomoge- 
neous corona is equal to that of a wedge corona at each annulus {R, R + dR). For spherical 
clouds of radius Rc{R), the coronal mass above (or below) the disk at each annulus is given 
by: 



2TTpoHRdR 



Tic-upcRtdR. 



which gives a clump radius 



Rc{R) 



R^ 



ndpc/po) 



1/3 



(9) 



(10) 



The covering fraction can be estimated by calculating the probability that a given region 
of the disk is not covered by a coronal cloud. Consider an annulus of area 2ttRAR. The 
expectation value for the number of clouds over it is = UcAR. If each cloud has cross 
sectional area irR^ ^ 27rR AR, this probability is given by 



P 



TlRl 



2'kRAR 



N 



N^oo 



exp 



2R 



exp 



2 " \2Rp, 



2/3' 



'11^ 



The covering fraction is simply fc = l—p. 

In Figure [8] we show a series of images of the BH accretion disk plus a clumpy corona, 
again sorted by photon history as in Figure [2j The corona is made up of 100 clouds dis- 
tributed randomly with scale height ratio H/R = 0.1 with constant probability per unit 
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radius (i.e., dN/dR = nc = 1) inside radius R = lOOM. The mean local optical depth 
To = 1, and the density contrast Pc/Po = 10. For the clumping algorithm described above, 
these parameters correspond to a covering fracti on in the inner d i sk of fr ~ 5%, comparable 



to that inferred for the hard state of Cyg X-1 by lPoutanen et al.l (119971 ). As in Figure [21 the 
BH has spin a/M = 0.9, mass lOM©, thermal flux Ltherm = O.lLEdd, corona temperature 
Tc = 100 keV, and observer inclination angle 75°. The direct flux from the disk, plotted in 
Figure [8^, clearly shows a number of distinct shadows where optically thick clouds block the 
disk. 

The reduced covering fraction of the corona leads to a greater flux from return radiation 
(Fig.[8)D), which otherwise behaves much the same as in the sandwich geometry, contributing 
a strong vertical component to the polarization signal. As in the direct image, the return 
radiation that eventually reaches the observer is also blocked in places by intervening coronal 
patches. In the lower panels of Figure [H we show the scattered photons (Fig.lHt: one scatter; 
Fig. [Hi: many scatters). Now the polarization of the hard flux is significantly less coherent 
than in the sandwich geometry. Many photons scatter multiple times within a single cloud, 
eventually emerging with no particular polarization direction. At the same time, the global 
geometry of the system is still essentially the same as the sandwich model, and photons are 
more likely to experience multiple scattering events by propagating from one coronal clump 
to another, constrained to move roughly parallel to the disk surface. As in the homogeneous 
model, this leads to a net vertical polarization in the observed flux at high energies. 

In Figure M we plot the flux, polarization degree, and polarization angle as a function 
of energy for the inhomogeneous corona and observer angles of 45°, 60°, and 75°. As in the 
sandwich geometry, here too we see that with increasing inclination, a greater fraction of the 
observed flux comes from scattered photons. Due to the smaller covering factor, a greater 
fraction of the thermal flux can escape at all inclinations, but the spectrum is still dominated 
by Comptonized photons above 7 — 8 keV. From equation (ITU]) , we see that the individual 
coronal clumps have typical optical depths of r 10 in the inner disk. This naturally leads 
to a harder power-law spectrum relative to that of the wedge geometry. 

The polarization signature for the inhomogeneous corona is similar to the sandwich ge- 
ometry at low energies, but noticeably different above a few keV. Below the thermal peak, a 
sandwich corona merely extends the classic Chandrasekhar atmosphere, leaving the polariza- 
tion unchanged. The scattering contribution (dot-dashed curves) is more weakly polarized 
when the gas is clumped because it is easier for a photon emitted at a large angle relative to 
the disk normal to scatter once off a coronal cloud and then reach an observeicl, contributing 



"^The clumpy corona can be thought of as a coUcction of hard scattering spheres where photons arc unable 
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a vertical component to the polarization, slightly diluting the dominant horizontal signal at 
low energy. On the other hand, the degree of polarization at high energies is significantly 
lower for a clumpy corona than for the sandwich corona. This is because many high-energy 
photons scatter multiple times within a single cloud, ultimately emerging with little or no 
net polarization due to the spherical symmetry of the cloud. 

Lastly, we note that when the corona is clumped, the shape of the transition from 
horizontal to vertical polarization more closely resembles that of the thermal state than that 
of a homogeneous sandwich corona: relatively smooth and in the ip < direction. Because 
of the small covering factor of the corona, thermal seeds emitted in the inner disk have a 
greater chance of returning to the disk or even scattering off a coronal clump on the far 
side of the BH, thus leading to a tp{E) distribution quite similar to that of the thermal 
disk with no corona (SK09). Not surprisingly, the polarization signature is sensitive to the 
compactness of the coronal hot spots, as shown in Figure [TOl In the limit of Pc/ Po = 1? 
the signal more closely resembles the homogeneous sandwich result: a sharp transition at 
2 — 3 keV, with ip > and the polarization degree rising above 6% at high energies (this 
case does not identically reproduce the homogeneous sandwich when is finite because, 
although the centers of the clouds must lie within H of the plane, their outer portions 
may extend farther, resulting in an inhomogeneous geometry even when Pc/Po = !)• As 
the compactness increases, the covering fraction and the overall symmetry decrease, giving 
a smoother transition and smaller amplitude above a few keV. In Figure [11] we show the 
dependence of the polarization signal on the number of clumps in the corona, holding the 
overdensity fixed at Pc/Po = 10. In the limit of A^'^ump oo, the homogeneous result is 
reproduced, and in the limit of A'^dump — we get the pure thermal result. 



SPHERE GEOMETRY 



The last model we consider is a spherical corona immediately surrounded by a trun- 
cated disk, a g e ometry motivated by a va ri ety of spectral and tim i ng observations (e.g . 



Gierlihski et al. 



(119971 ): iDone &: Zvckil fll999l ): lMakishima et all (120081 ): llngram et al.l (120091 ): 



^ — ^/ / ■ — — ui — -■ \i V ' ' — ■ \| ■/ 7 ' o — — ■ — ■/ 7 

Done fc Trigd (120091 )). The disk extends in to radius -Redge with a Novikov-Thorne emissivity 
profile, and the corona is defined by its temperature Tc and optical depth tq, as measured 
from the horizon out to Redge- The electron density in the corona is taken to be constant, and 
in keeping wit h its spherical shape , the fiuid is at rest in the frame of a zero-angular momen- 
tum observer (IBardeen et al.lll972l ). In order to most closely match the spectral properties of 



to forward-scatter, unlike the homogeneous atmosphere of diffuse electrons. 
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the sandwich corona, we also embed thermal seed photons with a uniform distribution inside 
the corona. The net flux and spectrum as a function of radius is determined by matching the 
shell-integrated flux to that of a standard Novikov-Thorne disk down to the ISCO. These 
coronal seed photons are emitted isotropically with zero polarization. 

In Figure [T2I we show the same photon-sorted images as in Figure [21 now for the spherical 
coronal model with -Redgo = 15M. The direct disk flux behaves as expected: because it is 
emitted from regions that are almost entirely in the non-relativistic regime, its polarization is 
very close to the Chandrasekhar value. There is now also a significant amount of unpolarized 
thermal flux coming directly from the coronal seeds. These seed photons can also scatter 
off the surrounding disk, in which case they are classified here as return radiation, again 
strongly polarized in the vertical direction. The corona-scattered flux shown in Figures 
[T2t,d closely resembles that of a scattering-dominated atmosphere around a star or planet: 
zero polarization from surfaces normal to the observer, and strongly-polarized (but limb- 
darkened) flux from the edges. However, the disk blocks a portion of the bottom half of 
the sphere, with the very bottom the most likely to be blocked. Consequently, the net 
polarization from a spherical corona is non-zero and oriented in the vertical direction. We 
have not shown here the contribution from photons that scatter in the corona and then off 
the disk, but these also are strongly polarized in the vertical direction, much like the image 
in the upper-right panel. 

Just as in the analogous figures for the sandwich and clumped models. Figure [12] shows 
that the total flux spectrum is dominated at low energies by the thermal contribution (here 
the sum of the disk proper and the thermal seeds embedded in the corona), while it is 
dominated at high energies by inverse Compton scattered photons. Because the coronal 
density and temperature were specifically chosen to give the same Compton y-parameter as 
in the sandwich geometry (see Fig. [3]), it is not surprising that the spectra agree so closely. 

The polarization signature is relatively straightforward to understand, especially in the 
context of Figure [12) the polarization from the direct radiation is almost exactly horizontal, 
decreasing steadily in amplitude with energy as a greater proportion of the observed flux 
comes from the unpolarized and higher-energy seeds in the corona. The scattered flux is con- 
sistently vertical, since even the lower-energy photons that have scattered only once have a 
net vertical orientation (see the lower- left panel of Fig. [T2|) . Because the unobscured coronal 
hemisphere looks much the same from each of the three viewing angles, the polarization am- 
plitude at high energies is only weakly dependent on observer inclination. Combining these 
effects gives the now-familiar behavior of horizontal polarization at low energies governed 
by the Chandrasekhar limit transitioning to vertical orientation above the thermal peak. 
The spherical geometry leads to a particularly sharp transition, at least for -Redge ^ lOM, 
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but it is still clear that a negative ip is preferred in the transition region. Again this can 
be understood by the geometry of the return radiation scattering off the region of the disk 
most strongly beamed towards the observer. In this case, the coronal photons originating 
from the central region above the disk plane replace the return radiation that similarly had 
passed near the photon orbit, just above the BH. 

In Figure [Ml we investigate the effects of varying Redge, while keeping the total optical 
depth of the corona constant. Since the thermal disk moves in with decreasing -Rcdgc? "we 
find the transition point moves to higher energies, as a greater fraction of the flux is coming 
directly from the disk and is therefore horizontally polarized. The shape of the transition 
grows broader with decreasing -Redge as relativistic effects become more important, rotating 
the polarization angle as in the thermal state. In the limit of large -Redge? where Newtonian 
physics dominates the problem, the system becomes scale-invariant and symmetry constrains 
the polarization at any given energy to be exactly horizontal or exactly vertical. In the limit 
of very small -Rcdgc, we reproduce the thermal result, plus a small power-law contribution to 
the spectrum at high energy due to the ultra-compact corona. 

Holding the corona radius constant at -Rcdgc = lOM, in Figure [T^ we show the effects of 
varying the optical depth and temperature of the corona. As in Figure [6] for the sandwich 
geometry, here too we keep the Compton ^/-parameter fixed so that the resulting spectra 
are nearly identical. However, unlike the sandwich model, for the spherical corona we find 
essentially no dependence on the coronal temperature or density. The reason is the essential 
simplicity of the scattering geometry: the polarization signature from the direct radiation 
alone is nearly the same, regardless of the coronal optical depth, and from Figure [13] we see 
that all scattered radiation is polarized exactly the same over the entire spectrum. Thus the 
behavior of the net signal is simply a function of the relative flux in each component, and 
this is held constant by fixing the Compton y-parameter. 



6. ACTIVE GALACTIC NUCLEI 



It has long been known that a significant fraction of the flux from AGN is emitted in 
the X-ray band (lElvis et al.lll978l ). As in stellar- mass BHs, this high-energy flux likely comes 
from lower-energy seed photons inverse Compton scattered in a corona of hot electrons with 
Tc ~ 100 keV. Similar to the stell ar-mass case, this leads to a relatively hard power-law 
spectrum with index a ~ 0.5 — 1 (iNandra et al.l Il99ll : iMushotzky et al.lll993l ). However, 
unlike the stellar-mass case, the temperature of the inner disk for an AGN will be well 
below a keV, leading to a thermal peak in the UV band. Furthermore, even when the 
disk is dominated by radiation pressure and electron scattering opacity, there should still 
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be a substantial fraction of metals that are not fully ionized, producing a large opacity for 
absorption above ~ 1 keV. 

Both of these AGN features — lower energy seed photons and an X-ray absorbing disk — 
lead to important differences in the polarization signature as compared to the stellar-mass 
models described above in Sections [3ll5l Because the seed photons start off with lower 
energies, they must scatter more times in the corona in order to reach the ~ 1 — 10 keV band. 
For a thin sandwich corona, this means that the scattering geometry is even more constrained 
than in the stellar-mass case, forcing the photons to move in a plane parallel to the disk 
surface, leading to a stronger vertical polarization. The AGN disk absorbs much of the 
incident X-ray flux from the corona, so the Compton ^/-parameter is effectively smaller than 
that of a stellar-mass system with the same coronal properties because scattering sequences 
are halted once a photon strikes the disk (we do not re-radiate the X-ray flux, but rather 
treat it as completely absorbed), thereby reducing the average path length the photons that 
escapes to infinity. An absorbing disk boundary condition with a sandwich corona therefore 
leads to an even higher degree of X-ray polarization because the photons are forced to scatter 
in a more constrained geometry before escaping the corona. 

In Figure [16] we show broad-band spectra from an AGN with central mass M = 10'^ Mq, 
spin parameter a/M = 0.9, thermal luminosity Ltherm = O.lLEdd, and observer inclination 
angle i = 45°. Four different corona models are considered: a homogeneous wedge with 
H/R = 0.1, To = 1, and = 100 keV, and three clumpy models with the same scale height, 
mean optical depth, and temperature, but with Uc = 1 and overdensities of Pc/po = 1, 
3, and 10. The inhomogeneous models can be characterized by their covering fractions as 
measured in the inner disk aX R = lOM; these overdensities correspond to fc = 0.25, 0.15, 
and 0.05, respectively. As can be seen in Figure [161 a greater covering fraction leads to 
a softer spectrum, while an inhomogeneous corona with smaller yet denser clumps gives a 
larger optical depth for those photons that do not escape directly from the disk, in turn 
giving a larger Compton y-parameter and harder spectrum. 

As mentioned above, the disk absorption in the AGN models results in a softer spectrum 
than stellar-mass BHs with the same coronal parameters. This can be seen by comparing 
Figure [TH] with the upper-left panels in Figures [3] (wedge geometry; fc = 1) and M (clumpy 
geometry; fc = 0.05). In both cases, the spectral index a is increased (i.e., softer) in the 
AGN case by about 0.5. The effect of this absorption on the X-ray polarization can be 
seen in Figure [T71 which plots the degree and angle of polarization as a function of energy 
for the same model parameters used in Figure [THl As expected, for the wedge geometry 
with fc = 1, the absorbing disk leads to a larger degree of polarization than the reflecting 
boundary condition used in stellar-mass BHs (~ 6 — 8% for AGN compared with the ~ 4% 
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shown in Fig. [3] for M = lOM© and i = 45°). 

For the clumpy coronae, however, we find the opposite effect: including absorption ac- 
tually reduces the degree of polarization. Recall from Section H] that inhomogeneous coronae 
generally produce a weaker polarization signal than the uniform- density sandwich geometry. 
This is because photons scattering multiple times in a single spherical cloud will ultimately 
escape with no net polarization. It is only through the global scattering geometry that a 
net vertical polarization is acquired, due to the small number of photons that scatter from 
one cloud to another, often reflecting off the disk surface at grazing incidence along the way. 
Since many of these photons are lost to absorption in the AGN case, global geometric effects 
play a smaller role, and thus the degree of polarization is diminished (~ 1% for AGN versus 
~ 2% for the stellar- mass case in Fig. Thus we flnd that for AGN, X-ray polarization is 
even more sensitive to the inhomogeneity of the corona than for stellar-mass BHs. 

From Figure [161 'we see that it should be possible to determine the AGN covering fraction 
from the X-ray spectrum alone, without polarization information. If so, then the degree of 
polarization could rather be used to constrain the inclination of the AGN disk, a parameter 
that can be quite difficult to constrain in many cases. In Figure (TH] we compare the degree 
of polarization for a range of inclination angles and two different corona models: wedge 
geometry with /c = 1 (solid curves) and a clumpy geometry with fc = 0.15 (dashed curves). 
Clearly the inclination will be easier to measure for a smooth corona, but even for a clumpy 
corona with fc = 0.15, a first-generation X-ray polarimeter should be able to distinguish 
between i = 15° and 45°. 

In the event that we can measure the covering fraction and inclination with other obser- 
vations, the polarization can constrain other corona parameters like the number of clumps 
Tie and their overdensity Pc/Po- From equation (ITT]) , we see that for a constant covering frac- 
tion, Uc oc {pcl PoYi so increasing the number of clouds both decreases their characteristic 
size and increases their density. Yet when averaged over a photon's entire path, a corona 
with larger Uc will appear more homogeneous, despite the increased density perturbations on 
very small scales. This in turn leads to a greater degree of polarization, as shown in Figure 
Uni Holding the inclination and covering fraction fixed at z = 45° and fc = 0.15, respectively, 
we vary the number density of coronal clumps. While all four cases plotted in Figure (TH] 
have nearly identical spectra, we see that polarization can distinguish them and thus give 
improved constraints on the corona geometry. The different models in Figure [19] would also 
likely have observably different timing properties: if individual clumps evolve coherently, but 
are independent of one another, systems with smaller Uc should produce greater amplitude 
flux variations in the X-ray band. 

In practice, AGN polarization measurements will have some additional observational 
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challenges not present in stellar-mass BH systems. First, the typical X-ray flux from nearby 
Seyfert galaxies and quasars is at least an magnitude smaller than the brightest galactic 
BHs. This means that much longer observation times will be required to reach a comparable 
level of polarization sensitivity. In addition, type 1 AGN are expected to have only modest 
inclinations (z < 45°) because at higher inclinations our view of the inner disk is likely 
blocked by the surrounding dusty torus (assuming the torus and disk are oriented in the same 
direction). There are also two potential sources of dilution. Radio- loud AGN generi c ally h ave 
somewhat larger ratios of X-ray flux to optical flux than radio-quiet IShen et al.l (120061 ). a 
fact plausibly interpreted as due to a jet contribution; in blazars, the jet contribution is 
almost certainly substantial. In addition, Fe Ka emission often accounts for a few percent 
of the 2-10 keV flux. Lastly, the much longer variability timescales for supermassive BHs 
make it unlikely that we could use the strategy discussed in Section [31 in which observations 
of multiple spectral states can jointly determine the orientation of the spin on the sky. 
This advantage would be absent in AGN systems, which for the most part do not undergo 
signiflcant state transitions on an observable timescale. 



7. DISCUSSION 

Using a Monte Carlo ray-tracing code in the Kerr metric, we have explored a variety 
of models for X-ray-emitting coronae attached to accreting stellar-mass black holes in the 
hard/SPL states or to radio-quiet AGN. Over a wide range of model parameters, the polar- 
ization swings from parallel to the disk plane ( "horizontal" ) at low energy to perpendicular 
( "vertical" ) at high energy. The location of this transition is generically at an energy a few 
times that of the highest temperature found in the disk; in the case of stellar-mass BHs, 

that means ~ 1 3 keV; in AGN, ~ 100 eV. The detailed properties of this transition, 

and the polarization amplitude at the low- and high-energy limits, provide the observer with 
information about the BH's mass, spin, inclination, and accretion rate, as well as the coronal 
properties: degree of (in)homogeneity, vertical scale height, temperature, optical depth, and 
covering fraction. 

We do not claim that a single or even multiple polarization observations of a given 
source will unambiguously measure all these model parameters. What does seem possible, 
even with relatively low precision data, is to rule out large regions of parameter space, and in 
so doing challenge a number of the traditional paradigms and toy models for BH accretion 
geometry. For example, if we were to observe an upper limit of 5 < 2% at ~ 10 keV 
from a stellar- mass system with known binary inclination i > 60°, we would have only two 
choices: to abandon any form of smooth corona geometry in favor of a highly asymmetric. 
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inhomogeneous accretion flow, or to posit a misalignment between the binary orbit and the 
inner accretion disk (thus giving lower inclination and polarization than expected). For 
AGN sources, where the expected disk inclination is lower [i < 45°), a null polarization 
measurement would rule out a smooth sandwich corona for anything but a nearly face-on 
system. 

As with any theoretical model based on a large number of parameters, we expect that 
certain physical properties of the BH system will be more easily constrained than others. 
For example, even if we assume prior knowledge of the observer inclination, BH mass and 
accretion rate, as well as the temperature, density, and geometry of the corona. Figure [5] 
suggests that, with a polarization sensitivity of 5 < 1%, one could only just distinguish 
between a Schwarzschild and extreme Kerr BH in the hard state. However, this also means 
that even if we have relatively poor constraints on the spin parameter (which is generally the 
case for most galactic BHs), we should still be able to measure other parameters robustly. 
By measuring the power-law slope of the spectrum around 10 keV, one could estimate the 
Compton y-parameter with reasonable accuracy, but still face a degeneracy between the 
corona temperature and optical depth. Again assuming a known inclination and corona 
geometry, we see from Figure E] that this degeneracy could be broken with a polarization 
observation. 

As another example, if the coronal properties are well-known (perhaps through more 
detailed spectral observations over a greater range of energies, which could allow an indepen- 
dent determination of the electron temperature), the polarization could be used to determine 
the disk inclination. Whi le the inclination could in principle be measured directly from the 



low-energy polarization (jConnors et al.l 119801 : iLi et al.l l2008l ). first-generation polarimeters 



are likely not going to be very sensitive below ~ 1 keV, where in any ca se, magnetic turbu- 



lence in the disk may reduce the net polarization via Faraday rotation (jPavis et al.ll2009l ). 
Yet as we see from Figure [31 a single high-energy polarization measurement at ~ 10 keV 
could give the inclination angle, assuming the other model parameters are known reason- 
ably well. In practice, it is more likely that we could use some independent method (e.g. 
optical light curves, Fe line fitting, etc.) to determine the disk inclination, and then use the 
amplitude of polarization in the ~ 1 — 10 keV band to determine the inhomogeneity of the 
corona. Coronae with large-amplitude density fluctuations produce weaker polarization at 
high energies and a more gradual transition from horizontal to vertical orientation, as shown 
in Figures [TO] and [TTl 



Poutanen et al.[ (119971 ) showed that, using basic physical arguments about the covering 
fraction of the corona and the relative flux in the hard and soft components, the smooth 
sandwich corona could be distinguished from a patchy corona on the basis of the spectral 
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properties alone. Yet in Section [5|, we have shown that the spectrum from a truncated disk 
around a spherical corona can be nearly indistinguishable from that of a wedge corona above 
an extended disk. However, the polarization signatures from these two models are quite 
distinguishable: the spherical geometry leads to a horizontal-vertical transition at lower 
energy and gives significantly weaker polarization at all energies. 

As seen in Sections [3] and HJ the sense of rotation in the polarization transition (i.e., 
i/j > or ip < 0) contains additional information that constrains the corona scale height and 
level of inhomogeneity. The angle of polarization at low energies defines the projection of 
the disk plane on the sky, but not the direction of the BH spin (as in many binary systems, 
the line of nodes is uniquely determined, but it is impossible to know which is the ascending 
node). From observations of the thermal state, we can use the sense of the rotation toward 
vertical polarization to determine the sense of orbital motion of the disk, thus defining the 
ip > direction. Given that result, the sign of ip in the coronal state can constrain the scale 
height of a wedge corona (Fig. [7]) or the covering fraction of a clumpy corona (Figs. [10] and 

m. 



Additional constraints on clumping properties pc/ Po and Uc can come from the amplitude 
of polarization above the thermal peak as follows: In stellar-mass systems, we might deter- 
mine the inclination of the disk from radial velocity observations of the c ompanion. In AGN, 



the in clination may be inferred from the blue edge of a broad iron line fiReynolds fc Nowak 



20031 ). In either stellar-mass systems or AGN, the slope of the hard power-law tail of the 
X-ray spectrum can give the covering fraction f^, which is related to Pc/po and Uc through 
equation (fTTj) . For a given i and /c, the degree of polarization at high energies can then be 
used to determine the coronal properties Pc/po and ric- 

Furthermore, with timing observations we expect an inverse correlation between the am- 
plitude of hard X-ray luminosity fluctuations and the number of coronal clumps, due simply 
to Poisson statistics. In a related way, polarization could provide valuable insight into the 
nature of poorly- understood timing phenomena such as quasi-periodic oscillations (QPOs). 
The QPOs observed in stellar-mass BHs are most prominent in th e hard X-rays (> 6 keV), 



even in the SPL state when most of the flux is at lower energies (IRemillard fc McClintock 



20061 ). suggesting that the QPOs come from the hot corona. By comparing the polarization 
signal with the timing properties at different epochs, we should be able to constrain the 
location, size, and coherence of the regions from which QPOs originate. For example, if we 
find that periods of large amplitude timing fluctuations correlate with low-amplitude polar- 
ization, then it is quite likely that both features are caused by a relatively small number of 
hot, dense clouds orbiting above the inner disk. 

Lastly, in Section [6l we showed how the special properties of AGN could be exploited 
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through X-ray spectropolarimetry. Because the seed photons in such systems are at much 
lower energies, more scattering events are required to inverse Compton scatter the seed 
photons to X-ray energies. Additionally, because the relatively cool disk is largely opaque 
to absorption below ~ 10 keV, the range of available scattering angles is more restricted. 
Together, these features particular to AGN systems lead to a stronger ability to discriminate 
between different coronal covering fractions, as well as the characteristic size and density 
of the scattering clouds. Combining spectral and polarization information, it should be 
possible to determine the coronal covering fraction and the disk inclination independently 
with moderate precision. In contrast to the galactic BHs, in the AGN case it is unlikely that 
we will be able to observe a single system in multiple accretion states, yet the much longer 
time scales should allow for more detailed comparisons between polarization and variability. 
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Fig. 1. — Schematic diagram of the disk and corona for the sandwich geometry. Thermal 
seeds are emitted from a thin disk in the midplane, then scatter off hot electrons in a corona 
with a wedge geometry and constant scale height H/R. The optical depth to electron 
scattering in the vertical is constant throughout the disk. For observers at high inclination 
(edge-on), the scattered photons will likely have high polarization, oriented perpendicular to 
the scattering plane, as indicated by the small black lines. Some photons emitted from the 
inner disk are deflected by the BH and then scatter off the disk on the far side, also leading 
to large amplitude polarization. 
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Fig. 2. — Ray-traced images of polarized flux from an accretion disk with a sandwich corona 
of scale height H/R = 0.1. The observer is located at an inclination of 75° relative to the 
rotation axis, with the gas on the left side of the disk moving towards the observer, which 
causes the characteristic increase in intensity due to relativistic beaming and boosting. The 
black hole has spin a/M = 0.9, mass M = IOMq. The observed intensity is color-coded 
on a logarithmic scale (normalized to the net intensity J du I^), and the energy-integrated 
polarization vectors are projected onto the image plane with lengths proportional to the 
local degree of polarization. The four panels correspond to the contributions to the observed 
flux (a) directly from the thermal disk; (b) return radiation scattered once off the disk; (c) 
photons scattered once in the corona; and (d) photons scattered multiple times in the corona. 
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Fig. 3. — Observed flux and polarization from an accretion disk with a sandwich corona 
geometry. The plots show the flux (left; arbitrary units of I'Fy), polarization degree (center), 
and polarization angle (right) as a function of observed energy, for inclinations of 45°, 60°, 
and 75° (top, center, bottom, respectively). The dotted lines represent contributions directly 
from the thermal disk, the dot-dashed curves are corona-scattered photons, and the solid 
curves are the total observed flux. 
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Fig. 4. — Degree and angle of polarization for a sandwich corona, varying the luminosity in 
the thermal flux. 
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Fig. 5. — Degree and angle of polarization for a sandwich corona, varying the spin of the 
BH, while holding fixed the total thermal flux. 
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Fig. 6. — Degree and angle of polarization for a sandwich corona, varying the optical depth 
and electron temperature, maintaining a roughly constant Compton-y parameter. 




Fig. 7. — Degree and angle of polarization for a sandwich corona, varying the scale height 
of the corona. 
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Fig. 8. — Same as Figure [2|, but for a hot spot corona with an overdensity factor of 10. There 
are roughly 100 hot spots distributed randomly within radius R < lOOM with scale height 
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Fig. 9. — Same as Figure [3], but for the hot spot corona geometry shown in Figure [HI 
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Fig. 10. — Degree and angle of polarization for a hot spot corona, varying the compactness 
of the coronal hot spots. 
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Fig. 11. — Degree and angle of polarization for a hot spot corona, varying the number density 
of coronal hot spots. 
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Fig. 12. — Same as Figure [2], but for a spherical corona surrounded by a truncated thermal 
disk with an inner edge at -Rcdge = 15M. 
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Fig. 13. — Same as Figure [3l but for the spherical corona geometry shown in Figure [T2l 
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Fig. 14. — Degree and angle of polarization for a spherical corona, varying the radius of the 
corona (inner edge of the disk) . 
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Fig. 15. — Degree and angle of polarization for a spherical corona, fixing -Redge = lOM and 
varying the optical depth and electron temperature. 
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Fig. 16. — Observed broad-band flux from a supermassive BH with M = lO^M©, a/M = 0.9, 
-Z^therm = O.lLEdd, and 1 = 45°. The different curves correspond to different covering fractions, 
as described in the text. 
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Fig. 17. — Degree and angle of polarization for a supermassive BH with the same parameters 
as in Figure [T6l 
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Fig. 18. — Degree of polarization for a supermassive BH with covering fractions fc = 0.15 
{dashed curves) and /c = 1 {solid curves), for a range of observer inclination angles. 
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Fig. 19. — Degree of polarization for a supermassive BH with covering fraction = 0.15, 
inclination i = 45°, and a range of clump sizes. From equation (fTTl) we see that, for a 
constant covering fraction, the number density of clouds ric is a function of the overdensity 
parameter Pc/Po- 
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